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Abstract: We discuss some of the implications of simulating QCD when the action used 
for the sea quarks is different from that used for the valence quarks. We present exploratory 
results for the hadron mass spectrum and pseudoscalar meson decay constants using im- 
proved staggered sea quarks and HYP-smeared overlap valence quarks. We propose a 
method for matching the valence quark mass to the sea quark mass and demonstrate it on 
UKQCD clover data in the simpler case where the sea and valence actions are the same. 
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1. Introduction 

To predict phenomenological quantities from lattice QCD with high precision requires the 
ability to simulate with light dynamical quarks. Ginsparg- Wilson fermions have the correct 
chiral and flavour symmetries. However, they are computationally expensive compared to 
improved staggered quarks. In the Nf = 2 + 1 improved staggered programme the square 
root of the fermion determinant is employed to reduce the number of dynamical flavours 
from four to two for the up and down quarks, and the fourth root is taken to reduce 
the number of flavours from four to one for the strange quark Q. Ensembles of gauge 
field configurations are then generated with these fractional power determinants as weight 
factors. There is no known local action to which this model corresponds. We define a mixed 
action as one where the action used to generate the ensemble of gauge configurations, or 
sea quark action, is different from the valence quark action used to determine hadronic 
observables on those configurations. Current Nf = 2 + 1 improved staggered simulations 
have a mixed action because the four-flavour staggered Dirac operator is used to generate 
the valence quark propagators rather than a local operator equivalent to that used in the 
ensemble weight. Unless a local operator can be found such that 

det Aocal = (det^ + m}) 1 / 2 (1.1) 

mixed actions are inevitable in the improved staggered programme. The Chebyshev poly- 
nomial approximation to the square root of (D st + m) is not the required operator as it 
has been shown to be non-local j|, |J]. That (D st + m) 1 / 2 is non-local does not imply 
that -Diocal does not exist, but serves as a warning, since the obvious candidate for such an 
operator fails. 

In the rest of this paper we assume that some -Di oca i exists so that the improved 
staggered ensembles are generated with an action in the same universality class as QCD. 
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We consider the case where the valence quark action is manifestly different from that 
of the sea and we choose the valence action that has the best chiral properties, that is, 
overlap valence quarks on an improved staggered sea. In || a local Symanzik action 
and the corresponding low-energy chiral effective Lagrangian are constructed for a general 
Ginsparg- Wilson valence action with Wilson sea quarks. Some of their considerations apply 
to more general mixed actions and, in particular, to overlap valence quarks on a staggered 
quark sea ||. 

Neuberger's overlap operator jj]] is given by 

DM = \ [1 + M + (1 - M) 7 5 e(#w(-p))] (1.2) 
where Hyj, is the Hermitian Wilson operator 

Hw{-P) = TsDw(-p) (1-3) 

with mass parameter < p < 2, and e(flw) is the matrix sign function of Hyj. The mass 
parameter fi is related to the bare quark mass am q through 

2p 

although we will ignore this below and write D ov (mo). The expectation value of some 
observable O in a model where the ensemble has been generated as 2 + 1 flavours of 
staggered quarks, with overlap valence quarks is 

(0) = ^J W(det{D st [C/]+m U(i }) 1 / 2 (det{D st [C/]+m s }) 1 / 4 e-^[ c/ ] (1.5) 

-Vi{ D o V [ U \( m i)}~ lr h 

where U are the gauge fields, Z is the partition function, {fji,r]i}, i = 1, • • ■ ,Nf, are the 
valence quark sources and S g is the gauge action. The real parts the eigenvalues of D ov 
are positive and bounded from below by the valence quark masses ntj, assuming rrii > 0. 
The expectation values are equal to those of a local field theory with action given by 

S = S g [U] + ^ XI (Aocal[tf] + m ud ) xi + Xs (Aocai[t/] + m B ) Xs (1-6) 

l=ud 

+ ^2{q l D m [U](m l )q i + D^U]^)^} 

i 

where the x fields are the one-component staggered sea quark fields, and the q fields are 
the overlap valence quark fields. The <j) fields are pseudofermion sea fields introduced to 
cancel the determinant of the overlap operator l| . For practical purposes the model can be 
regarded as having an exact SU(Nf\Nf)i,® SU(Nf\Nf)n<g>U(l)v symmetry when m{ = 
for i = 1, • • • , Nf Q. Restricting to transformations only in the valence quark sector, the 
infinitesimal chiral rotation is given by 

Sq = ierj 5 ^1 - TjAwJ <7 (1-7) 
6q = ieq ( 1 - ^D ov J 75 r 
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and possesses the correct U(1)a anomaly and an index theorem [10, [TlJ (for a review 
see pi). 

For Nf = 3, this model is in the same universality class as QCD when the sea and 
valence quark masses are matched. At non-zero lattice spacing, the separate chiral symme- 
tries for sea and valence quarks ensure that the lightest pseudoscalar meson mass vanishes 
at m va i = m sea = 0. This implies that the bare quark masses are related by 

m va i = C(a)m sea (1-8) 

where £ — > 1 as a — > 0. To date Nf = 2+1 simulations with staggered valence and fractional 
determinants of the staggered sea Q have set £(a) = 1. However, it is not obvious that 
this is the appropriate matching condition for overlap valence quarks on a staggered sea 
(or, for that matter, for staggered valence quarks). 

1.1 Matching the quark masses 

To match the sea and valence quark masses to their experimental values one would have 
to find an experimentally known hadronic state whose mass depends strongly on the sea 
quark mass. In principle, the if is one such hadron. The sea quark mass could be tuned 
until the rf has the correct experimental mass, whilst tuning the valence quark mass of the 
flavour non-singlet pseudoscalar meson to the pion. In practice, this is rather difficult, as 
the rf requires very high statistics calculations. An alternative would be to relate the bare 
sea and valence quark masses to each other via equation Jl.8|) , and then tune the flavour 
non-singlet mesons to their experimental values in the usual way. 

When the sea quark mass is infinite, i.e. quenched, then Bardeen et al. [IJJ] have 
demonstrated numerically that the model violates unitarity. This has also been observed 
numerically in partially quenched two flavour QCD and demonstrated in partially quenched 
chiral perturbation theory by fl5(| . We perform a similar analysis and show the same 
unitarity violation occurs when m va i < m sea . Our results suggest a criterion for matching 
the sea and valence quark masses. The quark masses can be tuned by varying the valence 
quark mass to see when these partially quenched pathologies appear for a given sea quark 
mass. This determines when the valence quark is lighter than the sea quark. 



Bardeen et al. 13, 14] show that the scalar correlator, 

Css(t) = e- fS mx)iP(x)^(0)ij(0)) (1.9) 

X 

is sensitive to this quenched pathology, because it couples to an rf — ir intermediate state. 
Shown in figure [l] are two of the diagrams which contribute to the rf propagator. Diagram 

a) , the "hairpin", has a negative coefficient. In full QCD, diagram b), with a series of vac- 
uum bubbles, cancels the effect of the hairpin diagram, so there is no negative contribution. 
In quenched QCD, only the hairpin diagram contributes, so the intermediate rf — tt state 
couples with a negative spectral weight. This gives the scalar correlator a negative value. 

In partially quenched QCD the situation is more complicated. The bubble in diagrams 

b) depends only on the sea quark mass, whereas the connected quark-flow lines depend only 
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Figure 1: Quark-flow diagrams contributing to the rj propagator. 

on the valence quark mass. Heuristically at least, the size of the contribution from diagrams 
b) can been thought of in the following way. When the sea quark mass is smaller than 
the valence quark mass, diagrams b) have a larger positive contribution than the negative 
contribution from diagram a). When the sea quark mass is heavier than the valence quark 
mass, diagram b) has a smaller contribution than a), which means the rf — n intermediate 
state couples to the scalar correlator with a negative weight. By monitoring the sign of the 
scalar correlator as the valence quark mass is varied it should be possible to match the sea 
and valence quark masses. 

To demonstrate this method, we examined the scalar correlator on the UKQCD Nf = 2 
clover data sets [16, 17], where the sea and valence quarks have the same action. So whether 
the valence quark mass is heavier or lighter than the sea is known. This data was generated 
with the Wilson plaquette gauge action and the clover quark action, where the coefficient of 
the Sheikholeslami-Wohlert term [O] was determined non-perturbatively fjl^| . For all data 
sets (3 = 5.2 and the volume is L 3 x T = 16 3 x 32. The values of the hopping parameter for 
the sea and valence quark masses, and the number of configurations are shown in Table |l|. 

The relatively poor signal-to-noise ratio for the scalar correlator implies the need for 
a large number of configurations. To improve the statistical resolution, we used a ratio 
of correlation functions, as the statistical fluctuations are correlated. In particular, we 
considered the ratio 

m - — c^j (L10) 

where PP denotes the pseudoscalar correlator. At sufficiently large times the ground states 
will dominate and R is then 

w- 1 - (L11) 

At the mid-point of the lattice 

R(T/2) = 1 - i^s e - Am3 7 2 (1.12) 
App 

where Am = ms - rap is the mass splitting between the scalar and pseudoscalar states. 
For a large enough lattice time extent, T, this ratio tends to unity at the mid-point. 

However, when the valence quark mass is lighter than the sea quark mass, the 7/ — tt 
state couples to the scalar correlator with a negative weight. Thus, a signal for the valence 
quark mass being lighter than the sea is R > 1. Figure |2| shows the ratio for different sea 
and valence quark masses. The open circle and filled square both have the sea and valence 
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quark masses equal, and R tends to unity at the mid point of the lattice. For the filled 
circles, R > 1 as t ~ T/2 at the 2a level, a signal for partial quenching, and indeed this 
data set has m va \ < m soa . This effect is clearly dependent on the sea quark mass, as the 
open and filled circles both have the same valence quark mass. 
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Figure 2: R(t) in equation (1.11) versus Euclidean time, t 
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Table 1: UKQCD dynamical clover (N f = 2) data sets for /3 = 5.2. 



Also shown in Table |] is the result of fitting equation ( l.llj) to the data. Clearly the 
ratio Ass/App is not a very well-determined quantity. However, it seems clear that this 
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ratio being negative is a signal at the 1 — 2a level that the data is partially quenched. 
Figure |3] shows both the scalar correlator and the ratio ( 1.1 1|) . At lighter quark mass 
and with fewer configurations, the fit results become rather dependent on the fit range 
chosen, but combining the fit information and examining these plots, it is clear that for 
ftsea = 0.13565, K va j = 0.13580 there is a signal for the negative weight state, and for 
fiVai = 0.13565 this signal is absent. 
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Figure 3: R(t) and Css(t) versus t for K sea = 0.13565. 

A precise matching of the sea and valence quark masses will be difficult to achieve, 
because the signal for the scalar ground state at large times for light quarks seems to 
disappear into the noise. When the valence quark mass is lighter than the sea, the signal 
for the negative weight rf — tt state is fairly strong. However, our results suggest that it is 
possible, in principle, (equivalently with very high statistics) to match the valence and sea 
quark masses. This is necessary to make sense of simulations with mixed actions when at 
least one of the sea or valence quark masses is out-with the chiral regime and matching to 
chiral perturbation theory is problematic. 

2. Overlap valence quarks on a staggered sea 



We have performed an exploratory study of overlap valence quarks on the MILC Nf = 
2 + 1 improved staggered configurations We measure the simplest states of the light 
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hadron spectrum, mesons and baryons, and the pseudoscalar decay matrix element for 
both light and heavy-light states. Due to a lack of computational resources, the number 
of configurations analysed was small. This prevented any realistic attempt at matching 
the sea and valence quark masses as described in the previous section. Whilst the results 
presented below, in figures 5 to 10, appear encouraging, the low statistics means they must 
be regarded as purely illustrative of the effectiveness of a mixed action approach. 

2.1 Smearing 

The overlap operator is only local for gauge configurations which are "smooth enough" ^l| . 
The MILC configurations we used have a lattice spacing of a ~ 0.125 fm and so are relatively 
coarse. Smoothing the gauge configurations should improve the localisation of the overlap 



operator. Moreover, smoothing the gauge fields by "HYP-smearing" [22] can improve 
the spectral properties of the Wilson-Dirac operator p3[| , which reduces the amount of 
computation required in the solver used to apply e(i/\y)- Indeed, HYP-smearing the gauge 
configuration does speed up the inversions. Furthermore, the low- lying eigenvalues of the 
staggered operator "mimic" the eigenvalue spectrum of the overlap operator when the 



configurations are smoothed in this way p3, 25, 26] suggesting that a smoothly behaved 



matching condition may exist for light quark masses. 

To examine the effect of multiple iterations of HYP-smearing, we studied the quark- 
antiquark potential on 624 quenched UKQCD configurations at f3 = 5.93 with a volume of 



16 3 x 32. The smearing parameters used were a± = 0.75, 02 = 0.60, and 03 = 0.35 [22]. 
Planar Wilson loops were used to extract the quark-antiquark potential, which was fitted 
to 

V(r) = V + ar - - (2.1) 
r 

Figure || shows the effect of multiple iterations of HYP-smearing on the string tension, a. 
Repeated HYP-smearing quickly altered the short-distance behaviour, while the medium- 
to-long distance behaviour remained relatively unchanged for a small number (< 3) of 
iterations. The effect on the potential of smoothing configurations has been studied many 



times before, following the work of Teper [27], and recently an extensive study for different 



actions and different smearings has been carried out [28]. Our limited study agrees with 
these previous results. We conclude for N < 3 the effect of smearing does not significantly 
alter the long range potential and thus the spectrum. 

2.2 The light hadron spectrum 

The overlap propagator calculations were performed on ten configurations from each of 
two ensembles produced by the MILC collaboration |2(J. One ensemble has am s = 0.05, 
ami = 0.03 and the other has am s = 0.05, ami = 0.02. Both have a lattice spacing 
a ~ 0.125 fm and linear size L ~ 2.5 fm. Three iterations of HYP-smearing were applied 
to each configuration. The overlap operator from the SZIN code [^] was then used to 
calculate propagators. These were created with seven different valence quark masses using 
the overlap multi-mass solver: four light and three heavy |3(J. Some of these results have 
been previously reported in [31]. 
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Figure 4: The effect of HYP-smearing on the long-range potential as measured by a. 



We performed simultaneous fits to three different correlators in order to extract the 
pseudoscalar meson mass (see figure |5|) . The fluctuations in the effective mass are larger 
than the apparent statistical errors, but this is probably due to underestimation of the 
variance on ten configurations. 

A partially quenched analysis was carried out, that is the sea quark mass was held 
fixed whilst varying the valence quark mass. Since we had multiple input valence masses, 
non-degenerate meson correlators could be constructed. Shown in figure ^ is the two- 
dimensional fit performed to (aMps) 2 versus valence masses m qi and m q2 , which allowed 
evaluation of the average u and d quark mass, m, from 

Ml = B(m qi +m q2 ) + A = 2Bm + A (2.2) 

where M n is the physical pion mass. The chiral symmetry of the operator should make 
(aMps) 2 vanish at zero quark mass. We do not constrain the fit to satisfy this condition, 
but within the limited statistics the parameter A is consistent with zero. This in turn 
allowed us to evaluate the strange quark mass, m s , from 

M\ = B (m s + rh) + A (2.3) 

where Mk is the physical kaon mass. 

We also determined the masses of the nucleon and delta baryon. The signal for the 
nucleon mass is very clean. Figure [?] shows the effective mass of the nucleon for the two 
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Figure 5: Pseudoscalar meson effective mass and simultaneous uncorrelated fit to three correlators 
[P = 975(7, = 97475?). The squares and diamonds are slightly offset horizontally for clarity. 

operators 

iVi(x) = Eijk^fC^)^ (2.4) 
N 2 (x) = e ijk (ipfC7 i >y 5 ip j )ip k . 

It is remarkable that we can see a signal for the negative parity partner of the nucleon on 
as few as ten configurations. This suggests that, despite their relative cost per propagator 
compared with staggered quarks, overlap valence quarks maybe the most cost effective way 
to extract precision light baryon physics from improved staggered configurations. 

Figure [§| shows the nucleon (upper plot) and decuplet (lower plot) masses versus the 
pseudoscalar meson mass squared. The lines are uncorrelated linear fits to the data. The 
values calculated by the MILC collaboration |2(| on their corresponding full ensembles 
are shown by open symbols. Both the nucleon and decuplet baryon masses from the 
overlap operator are significantly lower although, a priori, we don't know how to match 
the horizontal scales. The cut-off effects for the different formalisms will be different and, 
unless the matching function in equation ( |1.8| ) is very different from one, this suggests that 
the cut-off effects for the overlap baryons are be smaller. The nucleon mass shows some sea 
quark mass dependence, but the decuplet mass shows no variation. With ten configurations 
and relatively heavy sea quark masses, any trend is hard to spot. In the lower plot, the 
vertical dashed-dotted line shows the estimate of the rj sS mass squared, as measured by 
the overlap operator on these configurations. The horizontal dotted line is the physical £l~ 
mass in lattice units. Within large statistical uncertainties, this determination of the 
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mass at fixed lattice spacing agrees with the experimental value. Again, this may suggest 
that cut-off effects with overlap fermions are smaller than with staggered fermions, but, 
with data at only one lattice spacing, this remains speculation. 
The pseudoscalar decay constant, fps, is defined as 

Za^MPS) f0 
fps = M PS ■ (2 ' 5) 

We obtain Za from the axial Ward identity 

Z A {d^O) = 2m q {PO) (2.6) 
which we can express in terms of the pseudoscalar correlator, Cpp, and the pseudoscalar 



axial correlator, Cpa 4 - (Ol-AlPS 1 ) cancels in equation (2.5) and hence we require only Cpp 
to compute a renormalised fps- Once again, we performed a 2-d linear fit to the light 
non-degenerate pseudoscalars to calculate fps (see figure ||) and extracted the ratio of 
Ik/ fir (see table |2|). The value increases slightly with decreasing light sea quark mass in 
the right direction to agree with experiment. This is also evident from the slight change of 
the gradient with sea quark mass in figure ^. 

2.3 Charm Physics 

Heavy quark propagators essentially come for free in the overlap propagator calculation 
through the use of a multi-mass solver. However, lattice artefacts are 0(am q ) 2 and the 
heaviest input valence quark mass used is am q = 0.84, so (am,) 2 ~ 0.7. With the lattice 
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Figure 7: Nucleon effective mass for the heavier sea quarks, with etm q = 0.056 (equation |1 
square symbols show the negative parity excitation. 
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Sea Quarks 


fx/ fir 


Ids (MeV) 


om sea = 0.03/0.05 


1.03(3) 


226(14) 


am sea = 0.02/0.05 


1.08(4) 


232(11) 


Experiment |32j] 


1.22(1) 


266(32) 



Table 2: Pseudoscalar meson decay constants. 



spacing of a -1 ~ 1.5GeV, the calculation is at best on the limit of simulating charm. Due 
to the rapid decay in Euclidean time, we require double precision. However, this does 
not slow the solver down appreciably, as we need substantially fewer re-orthogonalisations 
against the projected eigenvectors of H\y in the linear solver than in single precision. 

These heavy quark propagators were used to calculate fo s (see table |||). The value 
of fos increases with decreasing light sea quark mass, in the direction of the experimental 



value, as can be seen from the change of gradients in figure 10 



The short distance behaviour of the potential has been altered by repeated smearing. 
As heavy quarks in quarkonium feel the short distance potential, this repeated smearing 
may be a source of worry. Indeed, examining the heavy-heavy correlator for the heaviest 
quark mass we do not see the effective mass reaching a plateau. It might be expected that 
a heavy-light state feels the effect of the short distance potential less. Indeed the effective 
mass for the heavy-light correlator reaches a plateau. This suggests that the heavy-light 
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Figure 8: The nucleon and decuplet baryon mass versus (aMps) 2 for the two ensembles. The open 
symbols show the baryon masses measured on the full ensemble with staggered valence quarks |20| . 
The lower plot shows the physical Q,~ mass in lattice units, where the lattice spacing is set by tq 
taken from the reference above. 

states are not suffering so much from the modified short-distance behaviour. In future 
work we anticipate using fewer iterations of smearing. 

3. Conclusions 

While staggered quarks offer the most cost effective way of simulating light dynamical 
quarks today, they require us to use a mixed action formulation of QCD. Outside the chiral 
regime of both valence and sea quarks, it is necessary to implement a matching procedure 
for the quark masses for the model to be in the same universality class as QCD. (Within the 
chiral regime, the partially quenched results may be matched to chiral perturbation theory 
and thence to QCD low energy constants.) Indeed, we show numerically that the partially 
quenched theory with m va i < m sea has similar negative metric pathologies to those observed 
by Bardeen et al. in quenched QCD. In principle, this observation provides a matching 
condition, but, just like the alternative approach of determining m sea by matching a flavour 
singlet quantity to experiment, suffers from poor signal-to-noise in practice. Despite these 
practical problems with matching, we obtain encouragingly good signals for flavour non- 
singlet hadron masses and decay constants using overlap valence quarks on a staggered sea 
quark ensemble. The potential gain from the simplicity of valence quarks with the correct 
flavour and chiral symmetries, together with the clean statistical signals, particularly for 



- 12 - 



0.25 



0.2 - 



> 

o 



0.15 




«3 
Oh 



0.1 



0.05 



• 


am = 

sea 


=0.02/0.05 




■ 


am = 

sea 


=0.03/0.05 




* 


kaon 






A 


pion 







_L 



_L 



0.1 



0.2 0.3 0.4 

(M p§ ) 2 (GeV 2 ) 



0.5 



0.6 



Figure 9: fps versus Mp S for the two ensembles. 

baryons, is good motivation for trying to improve on our exploratory attempts to match 
valence and sea quark masses. 
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